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Abstract -Most laser-induced femtosecond magnetism investigations are done in magnetic thin 
films. Nanostructured magnetic dots, with their reduced dimensionality, present new opportuni- 
ties for spin manipulation. Here we predict that if a magnetic dot has a dipole-forbidden transition 
between the lowest occupied molecular orbital (LUMO) and the highest unoccupied molecular or- 
bital (HOMO), but a dipole-allowed transition between LUMO+1 and HOMO, electromagnetically 
induced transparency can be used to prevent ultrafast laser-induced spin momentum reduction, 
or spin protection. This is realized through a strong dump pulse to funnel the population into 
LUMO+1. If the time delay between the pump and dump pulses is longer than 60 fs, a population 
inversion starts and spin switching is achieved. These predictions are detectable experimentally. 
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Introduction. — The discovery of femtosecond de- 
magnetization in ferromagnetic Ni pQ , which has attracted 
tremendous attention worldwide [2HZ], challenges tradi- 
tional wisdom and opens a new technological frontier for 
fast magnetic switching. The nonthermal inverse Faraday 
effect [8] is an excellent example, where one can manip- 
ulate the spin via the polarization of light, without in- 
volvement of the lattice. This has further inspired new 
theoretical and experimental investigations into the mech- 
anism of the femtosecond magnetism. Up to now, while 
the experimental focus has been on thin films and bulk 
materials [SJHU], magnetic nanostructurcs, with reduced 
dimensionality, unique symmetry and quantum confine- 
ment, are very attractive. For instance, their electronic 
and magnetic properties can be tailored systematically. 
Experimental studies in cobalt nanostructurcs show very 
pronounced spin excitations This presents a unique 
opportunity for electromagnetically induced transparency 
(EIT) in a spin system. EIT has been demonstrated in 
atoms and molecules |12j . and can be realized in some 
special designed quantum dots, with a proper state sym- 
metry. EIT in a three-level system relies on two coherent 
laser pulses: one as a control pulse and the other as a 
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probe pulse. The absorption of the probe depends on the 
control pulse (see fig. [1}. It would be fantastic if such 
EIT can be realized in a magnetic dot, since one can then 
fully integrate the light control into the magnetic storage 
on a femtosecond time scale, a daunting but technologi- 
cally significant task. This motivates us to pursue such a 
system. 

In this paper, we show that in a carefully designed mag- 
netic dot, EIT can be used to protect spin momentum re- 
duction. Using two laser pulses, one pump and one dump, 
we can not only protect a spin state but also switch it. The 
former is done through EIT, and the latter is achieved 
through population inversion. Spin protection does not 
occur naturally in any magnetic dots. The minimum con- 
dition for this to occur is that the optical transition be- 
tween the lowest unoccupied molecular orbital (LUMO) 
and the highest occupied molecular orbital (HOMO) must 
be dipole-forbidden, while the transition between HOMO 
and LUMO+1 or higher states is dipole-allowed. The spin- 
orbit coupling must not be zero. Our investigation is based 
on two different theoretical formalisms. We start with the 
first-principles calculation in ferromagnetic fee bulk nickel, 
and we then construct a model system with nearly identi- 
cal electronic and magnetic properties. The model system 
allows us to search the best candidate for spin protection 
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Fig. 1: Laser-induced spin protection and spin switching. A three-level system is shown. The transition between LUMO and 
HOMO is forbidden. 



Table 1: Eigenvalues, and spin and orbital momenta of fee 
Ni at fc-point of k = (101, 69, -7)/104. All the results are 
computed using the first-principles method. The Fermi energy 
is at eV. The fc-mesh is 104 x 104 x 104, so all the numbers 
are divided by 104. 
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£i (eV) 


S z (h) 


L z (h) 


1 


-4.065 


0.499 


-0.009 


2 


-3.740 


0.496 


-0.049 


3 


-3.691 


-0.494 


0.003 


4 


-3.112 


-0.498 


0.054 


■5 


-2.311 


0.496 


-0.024 


6 


-1.682 


-0.489 


0.019 


7 


-1.109 


0.490 


0.052 


8 


-0.565 


0.174 


-0.028 


9 


-0.499 


-0.174 


-0.001 


10 


0.351 


-0.499 


-0.017 


11 


4.863 


0.500 


0.002 


12 


5.086 


-0.500 


-0.003 


13 


9.873 


0.498 


-0.008 


14 


10.040 


-0.487 


0.007 


15 


10.261 


0.487 


0.009 


l(i 


10.434 


-0.498 


-0.008 


17 


12.152 


0.500 


0.001 


18 


12.252 


-0.500 


-0.001 



and switching. Future experiments can directly test our 
predictions. 

This paper is arranged as follows. In the second sec- 
tion, we present our theoretical scheme, followed by the 
laser-induced ultrafast demagnetization in the third sec- 
tion. The spin protection and switching is presented in 
the fourth section, and we conclude this paper in the fifth 
section. 

Theoretical formalism. — The first-principles 
method is certainly a method of choice for quantum dots, 
but there is a limit on what the method can do. It is 
difficult to isolate one physical quantity without affecting 
others in question. For instance, the spin-orbit coupling 
is a joint effect of the crystal potential and electron wave- 
function, and can not be parametrized easily at the first- 
principles level. To strike a balance between a real system 



Table 2: Eigenvalues (£;) and their spin (S z ) and orbital mo- 
menta (L z ). The spin-orbit coupling is chosen as A = 0.05 eV. 
The Fermi energy is set between ds and d§ in order to have 
small orbital and large spin momenta. 



i Orbital type £ t (eV) S z {h) L z (h) 
1 
2 
3 
4 
5 
6 
7 
8 
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10 
11 
12 
13 
14 
15 
16 



and a model system, we first perform a first-principles in- 
vestigation in fee nickel. While details have been presented 
elsewhere [T3], here in brief we use the full-potential aug- 
mented plane wave method within the density functional 
formalism as implemented in the WIEN2k code [2]. By 
solving the Kohn-Sham equation (in Ry atomic units), 

[-V 2 + V Ne + V ee + = EfMr), (1) 

we obtain the eigenfunctions and eigenenergies, from 
which we construct the spin and orbital matrix elements 
among all the band states. In cq. (TTJ), the first and second 
terms arc kinetic energy and Coulomb interaction between 
the electron and the nuclei, and V xc is the Coulomb and 
exchange interactions. The spin-orbit coupling is included 
via the second-order variational method. In table [TJ we 
show the eigenenergies, spin and orbital matrix elements 
at one crystal momentum point k. These elements are the 
basis for our model. 
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-2.496 


0.444 


-0.072 


(1-7 


-2.271 


-0.443 


0.073 


d 3 


-1.754 


0.479 


0.008 


<h 


-1.547 


-0.479 


-0.009 


ds 


-0.750 


0.497 


-0.003 


d e 


-0.550 


-0.497 


0.003 


in 


0.015 


0.350 


-0.013 




0.086 


-0.350 


0.013 


d 7 


0.248 


0.487 


-0.016 


ds 


0.452 


-0.486 


0.017 


dg 


0.973 


0.448 


0.068 


dio 


1.195 


-0.450 


-0.069 


P3 


1.440 


0.499 


-0.001 


P-i 


1.490 


-0.499 


0.001 


P5 


2.845 


0.357 


0.012 


Pa 


2.915 


-0.357 


-0.012 
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Fig. 2: (a) Energy spectrum. The results are obtained with A = 0.05 eV. Ef represents the Fermi energy. The pump laser is 
tuned to the transition between dg, and pi , while the dump to the transition between da and pi . (b) Ground-state spin angular 
momentum change as a function of A. (c) Ground-state orbital angular momentum change with A. (d) Transition matrix 
elements \fi a b\ between ds and pi and between ds and pi as a function of A. 



We introduce a sixteen-level model, with the spin-orbit 
coupling as an input parameter. The model consists of 
three p-orbitals (Yn, Flo and Yi_i) and five c?-orbitals 
(Y 2 2, Y 2 i, Y 20 , and F2-2) for each spin channel. 

Here the Y m 's are the spherical harmonics, and the radial 
part of the wavcfunction is included through the param- 
eters (see below). The overlap between p- and <i-orbitals 
is ignored to minimize the number of parameters used, 
but the overlap between d-orbitals and that between p- 
orbitals are introduced. To mimic a ferromagnetic ground 
state, these d-orbitals lie about 2 eV below the p-orbitals. 
Within these approximations, our Hamiltonian consists of 
two parts, Hq = Hd + H p . Here the Hamiltonian for d- 
states is 

See eq. [2] 

where Eda is the d-orbital energy with spin index a, I is the 
orbital angular momentum operator, and s is the spin an- 
gular momentum operator, a m(J (c mcr ) is the electron cre- 
ation (annihilation) operator, creating (annihilating) an 
electron in orbital m. We choose the spin-orbit coupling 
(SOC) constant A = 0.05 eV, except in figs. H», HJc) 
aiid^d). The Dirac bracket |2, m, a) in eq. © represents 
Y2m multiplied by the spin wavefunction \a). The last 
term, originating from the crystal and exchange fields, de- 



notes the overlap between orbitals with different to, with 
6 = 1 eV. (A similar Hamiltonian can be set up for the 
p-orbitals.) Based on our first-principles calculation for 
fee Ni [9JH3], we choose e<jt = —0.75 eV, e^j. = —0.55 
eV, e p -\ = 1.44 eV and e P i = 1.49 eV. While the precise 
orbital energies are not critical here, in order to realize a 
strong spin momentum, we find it is necessary to adjust 
the overlap term if different orbital energies arc used. An 
extension to /-orbitals is straightforward. 

We first compute the spin, orbital and dipole transi- 
tion matrices by diagonalizing the Hamiltonian. The or- 
bital characters of the eigenstates are determined by their 
eigenvectors. The eigenvalues, spin and orbital momenta 
are tabulated in table [5] The six lowest eigenstates are 
all d-states, followed by two p-states. It is clear from the 
table that to realize large spin but small orbital momenta, 
the number of occupied states is not arbitrary. We choose 
the Fermi energy between d§ and d^ so the five lowest 
(i-orbitals are occupied. The resultant ground-state spin 
momentum is (s z ) = 0.498?i, and the orbital momentum 
is (l z ) = — 0.0027?i. This result can be compared with our 
first-principles results in table [TJ and we find that they 
match very well. For instance, the majority of spin values 
are around 0.5?i and the orbital momentum is below O.lfi.. 
This gives us confidence that our model is able to simulate 
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Time, t (fs) Time delay, T (fs) 

Fig. 3: (a) Population change as a function of time upon laser excitation. Only the three states pi, P2 and d$ are shown (the 
population in all the other states is nearly unchanged). Inset: Laser pulse, (b) Spin momentum change as function of time. 
The oscillation is due to the interference between states p\ and P2- (c) Orbital angular momentum change, (d) Population 
change excited by pump and dump pulses as a function of the time delay between the two pulses. Electromagnetically induced 
transparency (EIT) occurs for the transition between d$ and pi when the delay is close to fs, and population inversion occurs 
for the delay longer than 60 fs. (e) Change from spin protection to spin switching, where the time delay between the pump and 
dump pulses is longer than 60 fs. 



a ferromagnetic system. 

Figure [UJa) shows the energy spectrum with the rele- 
vant transition states labeled. Figures [2jb) and[2jc) show 
the dependence of the spin (s z ) and orbital angular mo- 
menta (l z ) on the SOC constant A. It is interesting to 
note that the spin decreases with A, but we also find this 
to be parameter dependent, and in some cases the depen- 
dence is not monotonic. The orbital momentum is zero 
if A is zero, or complete quenching, but with nonzero A 
the quenching is not complete. This reproduces the well- 
known fact that the SOC drags some orbital momentum 
with it. We find transition matrix elements are larger 
between two states with a similar spin moment, which 
independently validates our prior first-principles results 
P21[T5]. For instance, the transition matrix element (z- 
direction) between d 5 (s z = 0.497?i) and p\ (s z = 0.350?i) 
is 0.196 £, while the element between d 5 (s z = 0.497?i) 
and P2 (s z = — 0.35QK) is -0.078 £, where £ is a constant 



resulting from the radial part of the wavefunction. This 
demonstrates again that our model is applicable to ferro- 
magnetic systems. Therefore, even with SOC, the optical 
transition still prefers spin-conserved transitions over spin- 
nonconserved transitions, though for the demagnetization 
process, the latter transitions are most important. The 
systematic change of the transition matrix elements with 
the SOC for the above two transitions is shown in fig. 
[21(d) . To our knowledge, this is the first numerical result 
that shows clearly the transition matrix element drops for 
a "spin-conservation" transition (here between d§ andpi), 
while the elements for the spin-flip transition (between d$ 
and pi) increase. This demonstrates the critical role of 
the SOC in ultrafast demagnetization. 

Laser-induced coherent magnetization. To 

simulate the ultrafast spin evolution, we introduce a laser 
pulse with Ei{t) = Ai cos(wii) exp[— (t — ti) 2 /rf], where 
Ai, Ui, ti and Tj are the laser amplitude, frequency, time 
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delay and pulse duration of pulse i, respectively, and t is 
time. The laser has a Gaussian shape, with duration of 40 
fs (see the inset of fig. Ela))- We tune the laser frequency 
to be resonant with the dipole-allowed transition between 
g?5 and p± . The dynamic simulation starts with solving the 
Liouville equation for the density matrix p as [41 I131IT5] . 

ih^ = [H,p], (3) 

where H consists of the Hamiltonians Hd and H p for the 
d- and p-states and the interaction Hi between the sys- 
tem and laser field. Hi = —eE(t)J2ij(i\ z \3) c l c ji where 
— e is the electron charge, the light polarization is along 
the z-direction, and is the transition matrix element 

between states i and j. The merit of our approach should 
not be underestimated. Although the Liouville equation 
is equivalent to the time-dependent Schrodinger equation, 
the Liouville equation naturally takes into account the 
Pauli exclusion principle and includes the antisymmetry 
of the many-body wavefunction. However, most of the 
time-dependent Schrodinger equations from textbooks are 
for a single particle only, which is not suitable even for a 
non-interacting case. 

The general form of the density matrix is a convoluted 
integration equation. We express the density matrix in a 
state representation as 

ifr = ^" _ £rn)Pnm + ^(-^nfePfcrn - PnkH[ m ), (4) 

k 

where k, n and m are the state indices, £ n is the eigenen- 
ergy of state n, p nrn is the density matrix between states 
n and m, and H 1 is the interaction between the laser field 
and the system. If we integrate over time, we have 

Sec cq. [5] 

which is the general form of the density matrix. It is 
clear that the time evolution of the density depends on 
its history and the laser field. Further simplifications are 
possible, only if we know the profile of the laser field. For 
instance, if the laser is a continuous wave and is very week, 
the density matrix can be approximated by the first or- 
der density matrix whose time evolution is proportional 
to exp[— i(uJiaser — w nTO )t], where uj nm is the transition 
frequency between states n and m, u>i aser is the laser fre- 
quency, and t is the time. If the laser is a pulse, then the 
profile of the density matrix sensitively depends on the 
laser pulse duration, laser frequency and the electronic 
states involved as can be seen from the above equation. 
In this case, an analytic solution is generally not possible, 
and a numerical integration is a must. 



Figure EJa) shows that the population change closely 
follows the profile of the laser field. Upon the laser ex- 
citation, the c^-state loses its population from 1 to 0.25 
within the first 40 fs. In the meantime, the pi-state gains 
the same amount of population from to 0.75. Other 
states have a very small change. For instance, p2 only has 
a tiny increase around fs. Once the laser field is over, the 
population or diagonal element of the density matrix pa 
change is stabilized. However, this does not mean that all 
the elements of the density matrix become time indepen- 
dent. The interference between states, which is induced 
by the laser field initially, remains. This is precisely what 
happens with the spin momentum change. 

Figure [3fb) shows the spin change as a function of time. 
The spin drops similarly within the first 40 fs, followed 
by an oscillation with period 57.7 fs. The first drop re- 
sults from the spin momentum difference between d§ and 
Pi, as seen from table [5J where (d 5 \s z \d 5 ) — 0.497?i and 
(Pi\ s z\pi) = 0.350?i. The oscillation is from the interfer- 
ence between p\ and p 2 , which can be verified by switching 
on/off the density matrix element p PljP2 - It is remarkable 
that even though the population of p 2 is almost zero, its 
interference with p\ is directly responsible for this salient 
change. The period 57.7 fs corresponds to 0.071 eV, which 
matches the energy gap between p\ and p 2 exactly. The 
orbital momentum change is shown in fig. EIc). In con- 
trast to the spin change, the orbital angular momentum 
increases (more negative). Besides its initial drop, the or- 
bital momentum also has smaller rapid beatings. 

Spin protection and switching. As mentioned 
in the Introduction, the key to the electromagnetically in- 
duced transparency relies on the state symmetry. In the 
following, we demonstrate that one can not only reduce 
and flip a spin (spin switching), but also protect it (spin 
protection). We employ two laser pulses, one pump with 
field E\(t) and one dump with E 2 (t) [T7]. The laser pa- 
rameters of the pump are the same as above. The dump is 
tuned to the transition between p\ and de (see fig. EJa)), 
with duration of 60 fs and field amplitude of 2.5 times 
the pump. These parameters are carefully chosen to max- 
imize the effect and are fixed below. The only variable 
is the time delay T = t\ — ti between the pump and 
dump. Figure 0(cl) shows the populations as a function 
of T for three states, with the solid line for the originally 
occupied ds-state, dashed line for the originally unoccu- 
pied pi-state, and circle line for the originally unoccupied 
rig-state. Negative time delay means the dump proceeds 
earlier than the pump. The results in figs. [3jd) and 3(e) 
are computed after t = 200 fs. If the results oscillate with 
time, we show the time-averaged data. At T = —60 fs, 
similar to the above single-pump-pulse excitation, d§ loses 
population to pi, but now the dg-state gains some popu- 
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lation, in sharp contrast to the above single-pulse excita- 
tion. When T decreases (less negative), the population in 

increases sharply, but the populations in p\ and d$ drop 
sharply. In other words, the system becomes more difficult 
to excite and much less absorptive, or transparent to the 
light field. This is a manifestation of electromagnetically 
induced transparency (EIT) |18j . though in EIT normally 
cw lasers are used [TS]. What is novel here is that a time 
delay can induce EIT, which has a direct consequence in 
the spin change. Figure [3jo) illustrates that during EIT, 
the spin is mostly unchanged, i.e., spin protection. What 
is even more interesting is that when we increase the time 
delay over 60 fs, population inversion starts. 

Going back to fig. Eld), we see that the population in de 
becomes largest and that in dr, drops below both d§ and p\ , 
a total population inversion. Now if we monitor the spin 
change at the same delay, we find that s z changes its sign 
from positive to negative (see fig. [3je)), or spin switching. 
We expect our predictions are detectable experimentally 
in magnetic systems whose electronic and magnetic struc- 
tures are similar to those of our model. The key to the 
successful realization of EIT is that the dipole-transition 
between the LUMO and HOMO must be forbidden and 
the transition between HOMO and LUMO+1 or higher 
states is allowed. Very recently, Hiibner et al. demon- 
strated the optical spin manipulation for magnetic logic 
operations in the Ni3Na2 cluster P21I20]. Therefore, the 
possibility to observe our prediction is very high. 

Conclusions. — We have demonstrated a new elec- 
tromagnetically induced spin protection and switching in 
a specially designed magnetic dot. The minimum require- 
ment for such a dot is that the optical transition is al- 
lowed between the HOMO and LUMO+1, but not HOMO 
and LUMO. A strong dump pulse couples LUMO and 
LUMO+1, and directly controls the spin dynamics. As 
a result, the spin becomes harder to excite, i.e., spin pro- 
tection. We have also shown that if the time delay be- 
tween the pump and dump laser pulses exceeds 60 fs, a 
spin switch from spin up to spin down occurs. We expect 
both these predictions are detectable experimentally. 
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